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Abstract The accelerated Kepler problem is obtained by adding a con- 
i-C ■ stant acceleration to the classical two-body Kepler problem. This setting 

models the dynamics of a jet-sustaining accretion disk and its content of 
forming planets as the disk loses linear momentum through the asymmet- 
ric jet-counterjet system it powers. The dynamics of the accelerated Kepler 
problem is analyzed using physical as well as parabolic coordinates. The lat- 
ter naturally separate the problem's Hamiltonian into two unidimensional 
Hamiltonians. In particular, we identify the origin of the secular resonance 
in the accelerated Kepler problem and determine analytically the radius of 
^ ' stability boundary of initially circular orbits that are of particular interest to 

the problem of radial migration in binary systems as well as to the truncation 
of accretion disks through stellar net acceleration. 

(N ■ 

Keywords Kepler problem • stellar jets • extrasolar planets • stellar binaries 

o . 

t"*** , 1 Introduction 

o 

The recent studies of the consequences of the acceleration induced by stellar 
jets on extrasolar planets and protoplanetary disks (Namouni 2005, 2007, 
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hereafter Papers I, II; Namouni and Zhou 2006) have shown that the basic 
dynamics of these systems rely on the accelerated Kepler problem (hereafter 
AKP). Although it is uncommon in celestial mechanics, the problem obtained 
by adding a constant acceleration to one of the two bodies in an r _1 poten- 
tial is at the basis of the Stark effect in quantum mechanics (Brandsen and 
Joachain 1989). In this case, the acceleration models the constant electric 
field applied to an atom. 

An interesting property of the AKP is its integrability. This result was ob- 
tained by Epstein (1916) who used the Hamilton- Jacobi appraoch in terms of 
parabolic coordinates to separate the AKP into two unidimcnsional Hamil- 
tonians (Landau and Lifschitz 1969). The integrability of the AKP in the 
position-velocity space was expressed by Redmond (1964) who derived a 
generalized Runge-Lenz invariant for the Stark effect. This invariant together 
with the energy and the angular momentum component along the direction 
of acceleration are the three integrals of the AKP. Recently, it was recog- 
nized that some of the dynamical consequences of time-variable jet-induced 
acceleration on the evolution of accretion disks and the forming planets they 
contain can be ascertained from the AKP (i.e. from the dynamics associated 
with a time- independent acceleration). To this end, analytical derivations 
in terms of the physical variables (positions, velocities, orbital elements) as 
well as perturbation theory through secular averaging and sudden impulse 
methods were used to determine the eccentricity and inclination excitation 
of extrasolar planets (Paper I) as well as the vertical profile of jet-sustaining 
accretion disks (Paper II) . 

The use of parabolic coordinates has the advantage of reducing the AKP 
into two integrable unidimensional problems. This however comes at a price: 
parabolic coordinates are not physical. Consequently, the aim of this paper 
is to link the two analytical approaches of perturbation theory and parabolic 
coordinates in order to give a canonical interpretation of the physical prop- 
erties described in the works mentioned above. 

The paper is organized as follows: in Section 2, we show how the AKP 
arises in Keplerian systems with asymmetric momentum loss. In Section 3, 
we compare some of the AKP's simple qualitative features to those of the 
classical Kepler problem. Section 4 is devoted to the AKP's formulation in 
parabolic coordinates. The secular resonance at the basis of the eccentricity 
excitation of extrasolar planets is studied in section 5 and its origin in the 
non-averaged problem is given in terms of parabolic coordinates. In section 
6, we derive the dependence of the Keplerian boundary radius on the initial 
inclination for circular orbits. Section 7 contains further directions of research 
in the AKP. 



2 Prom the classical Kepler problem to the AKP 



Consider a two-body Kepler problem where the larger mass object is loosing 
linear momentum at the constant rate Pi OS s ■ The equations of motion of the 
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two bodies in an incrtial frame are given as: 

dv GM(x-X) dV Gm(x-X) • 

where m, M, x, X, v, V are the masses, positions and velocities of the two 
bodies and G is the gravitational constant. In terms of the center of mass 
and relative motion coordinates, x r = x — X, (to + M)x g = rax + MX and 
the corresponding velocities, v r , v g , the equations of motion are written as: 



dv 9 Pi oss dv r G(m + M) 



dt (to + M) ' dt |x r | 3 M 



(2) 



These equations show that the motion of the center of mass is accelerated 
in the inertial frame by an amount equal to the ratio of the momentum loss 
rate to the total mass of the system. The relative motion of the small body 
remains independent of the center of mass's motion while it is accelerated in 
the opposite direction to momentum loss by the ratio of the momentum loss 
rate to the central mass. 

Calling A = — Pi oss /M the constant acceleration of the relative motion 
and dropping the index r for simplicity, the relative motion of the smaller 
object around the larger one is given by: 

dv k 

■x + A, (3) 

where k = G(M + m). This equation defines the accelerated Kepler problem. 
It describes the dynamical evolution of a planet embedded in a jet-sustaining 
disk that looses mass and linear momentum from its inner parts (the so-called 
jet launching region, Paper II). 



3 Constants of motion and equilibrium points 

As the perturbing acceleration is constant, the dynamical problem described 
by © remains conservative and amounts to the addition of the potential 
R = A • x to the two-body Kepler problem. This in turn gives the energy 
integral of the AKP as: 

The acceleration A defines a preferred direction in the AKP. As a result, the 
projection of the specific angular momentum h = x x v along this direction, 

h t = h- A/ A, (5) 

is conserved as h = x x A. The third integral (Redmond 1964) is given as: 



(6) 
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Fig. 1 Zero- velocity curves of the classical Kepler problem in the cylindrical coor- 
dinate plane (p,z). The unit distance is the Astronomical Unit (AU). The central 
object's mass is equal to the Sun's and the value of h z used in the expression of W 
is that of a planar circular orbit at p = 36 AU. 



where e = v x h/k — x/|x| is the Runge-Lenz vector. The constant (3 replaces 
the Runge-Lenz vector of the classical Kepler problem. 

To further compare the AKP with the classical Kepler problem, we exam- 
ine its orbits of least energy. As the acceleration imposes a preferred direction, 
we choose a reference frame where the z-axis is directed along A. Using the 
conservation of h z , the energy integral (j4|) can be written as: 

E=£±^ + W(p,z), where W(p,z) = ^- * -Ax, (7) 

and where p and z are the usual cylindrical coordinates in the plane orthog- 
onal to A. For a given value of h z , the zero- velocity curves (level curves of 
W) are shown in Fig. 1 for the classical Kepler problem and in Fig. 2 for the 
AKP. For the classical Kepler problem, only a stable equilibrium associated 
with circular motion is present. In the case of the AKP, the addition of the 
acceleration potential and its odd-symmetry introduce a new unstable equi- 
librium as well as a motion separatrix. The two equilibrium orbits (stable 
and unstable) of the AKP are located above the reference plane, z — 0, that 
contains the central object. The separatrix denotes the boundary where the 
acceleration from the central object matches the amplitude of the perturb- 
ing acceleration. Motion outside the separatrix is unbounded. To determine 
the location of the equilibrium points, we minimize the energy integral with 
respect to p, z, p, and z. The first two derivatives yield p — z — showing 
that least energy orbits are circular and planar. The derivative with respect 
to p yields the orbital rotation rate as a function of the projection of the 
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Fig. 2 Zero- velocity curves of the AKP in the cylindrical coordinate plane (p, z). 
The unit distance is the Astronomical Unit (AU). A constant vertical acceleration 
of amplitude A — 2 x 10 -10 kms -2 was added to the classical Kepler problem of 
Fig. 1. The odd-symmetry of the acceleration's potential removes the reflection 
symmetry of the Kepler problem with respect to z. The region of bounded motion 
is restricted by the appearance of the unstable orbit at (p, z) = (102.8, 103.5). 



angular momentum on the direction of acceleration, h z . Using the fact that 
orbits are planar and circular, h z = np 2 and the orbital rotation rate: 



k 



(p 2 + z 2 ) 



2\3/2 



(8) 



The derivative with respect to z amounts to equating the vertical accelera- 
tion, A, and the vertical pull of the central object and yields 



k z 
T 



(9) 



The circular orbits hover above the central object increasingly higher with 
distance. In other words, a circular orbit in the AKP is still planar but the 
corresponding orbital plane no longer contains the central object. The latter 
becomes increasingly distant from that plane as the orbital radius increases. 
Figure (3) shows the family of circular orbits of the AKP in the pz-plane. 
The uppermost orbit is located at p = and z — y/k/A. 

The locus of least energy orbits given by equations ([8|) and ([9|) does not 
indicate where the stable orbits end and the unstable orbits start. To find 
this limit, we express the specific angular momentum, h z , for the least energy 
orbits |9|) as a function of z. As the equilibria of W in the pz-plane are 
determined by the same h z , their coordinates are found by solving for z at 
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Fig. 3 Least energy orbits of the AKP in the pz-plane. The system's parameters 
are those of Fig. 2. The solid line stands for the stable equilibrium orbits (the 
sombrero profile) and the dashed line for the unstable equilibrium orbits (see Fig. 
4 for more details). 



constant h z and using equation Q to determine p. The maximum of h z {z) 
is reached at: 



i i 

(10) 



This location determines the last stable planar circular orbit of the AKP. 
Figure (4) shows the function h 2 z (z) as well as p(z) ([9]). The maximum of h z 
is where the stable and unstable points meet; above this value of h z motion 
is unbounded in all space. For z < z cr j t , planar circular orbits are stable 
and their profile in the pz-plane has the shape of a sombrero curved along 
the direction of acceleration (Fig. 3). In Paper II, the sombrero profile was 
identified as the instantaneous state of least energy of a jet-sustaining disk. 
This profile was used to ascertain the changes in the temperature, hydrostatic 
equilibrium and heating of the disk. 

Lastly, we remark that the existence of the constant of motion (3 shows 
that the separatrix of W is pendulum-like and does not involve any kind 
of chaotic motion. Orbits that are set up near the unstable point such as 
that in Figures 7 and 8 of Paper I tend to escape after a long time (relative 
to the orbital period) only because of the discrete nature of the numerical 
scheme used to integrate the orbit. In section 6, we discuss further the outer 
boundary of initially circular orbits whose orbital plane is inclined to the 
direction of acceleration such as those studied in Paper I. Before doing so, 
we first introduce the formulation of the AKP in parabolic coordinates that 
will be needed in sections 5 and 6. 




27A 



and 



Pcrit 



27A 
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Fig. 4 The identification of the equilibrium orbits' stability through the vertical 
component of angular momentum. The system's parameters are those of Fig. 3. 
Least energy orbits (lower maximum curve) are shown in the zp-plane (left scale). 
The function h 2 (z) (higher maximum curve, right scale) is determined from the 
equilibrium orbits (Eqns. [8] and [9} and normalized to the angular momentum of a 
circular orbit at 1 AU. The function h 2 z (z) yields two values of z for a given value 
of h z - These are the stable (solid line) and unstable (dashed line) equilibria of Fig. 
2. The maximum of h 2 (z) denotes the outer edge of the sombrero profile. 



4 The AKP in parabolic coordinates 

Parabolic coordinates (£, rj, 9) are related to the cylindrical coordinates (p, z, 9) 
through (Landau and Lifschitz 1969): 

*=\{£-ri), P=Vtv- (11) 

In terms of the orbital radius r 2 = p 2 + z 2 , these coordinates are written as: 

£ = r + z, rj = r — z. (12) 
The conjugate momenta of parabolic coordinates are given by: 

PZ = ^(Z + V), Pv = ^(t + V), Pe = (,r 1 6. (13) 

The Hamiltonian of the AKP can now be written in terms of the parabolic 
coordinates and momenta as: 

f = 2 M±^ + j!_^_a (14) 

where we used the relations ((TTj) and (fL2| in the expressions of the gravi- 
tational potential as well as the acceleration potential R = Az (third and 
fourth terms on the right hand side) . As the Hamiltonian is autonomous and 
independent of 9, we recover the two constants of motion: the energy |4|. 
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H = E, and the angular momentum component along the direction of ac- 
celeration ([5]) , po = h z . By applying the Hamilton- Jacobi approach with the 
separation of variables method, it can be shown (Epstein 1916, Landau and 
Lifschitz 1969) that the Hamiltonian (|14| decouples into two unidimensional 
Hamiltonians that are related by the energy E, the vertical component of 
angular momentum, h z and a third integral, (3 defined in equation (O, and 
given by: 

h 2 An 2 h 2 AF 2 

P=2 V p 2 r , + ^ + ^--Er 1 -k, -(3 = 2^l + ^-^--E^-k. (15) 

The two unidimensional Hamiltonians are therefore: 

9 h 2 A ri (3 + k s 
E = 2p ' + 2t 2+ ^-— ^ 

E = M + §-§ + ^- (17) 

The choice of vertical direction for the coordinate frame along the direction of 
acceleration implicitly imply that A is positive. Close examination of (I16p and 
(fl7| shows that it is (fT6| that has a potential function (last three terms on 
the right hand side) that has a minimum regardless of any constants because 
of the positive linear term associated with A. When motion is possible, orbits 
are bounded. Motion is possible only if E is larger than the potential function 
minimum corresponding to a value that wc shall call rj s . The motion in the 
r^-plane is shown in Figure (5). 

Motion in the £j>£-plane depends on the constants of motion as the neg- 
ative linear acceleration term needs to be compared to the other potential 
terms that appear in the potential function: 
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Fig. 6 Contour plots of the Hamiltonian (|17|l in the £pj-plane with the parameters 
h z = 1, A = 0.1, 0=1 and k — 2. The separatrix passing through the unstable 
point at £ u divides space into bounded and unbounded motion. 



The derivative of shows that the extrema satisfy the polynomial equation: 

+2 (/J- k)£ + 2h 2 z = 0. (19) 

For f3 > k, this equation has no solution (as A is positive) and Wj is a 
monotonic decreasing function. Consequently, motion is unbounded. For /3 < 
k, there are extremum points only if: 



hi < 



2 ( k ~ P) 



A~*. 



(20) 



In this case W has a relative (stable) minimum and a relative (unstable) 
maximum that we call £ s , and £ u respectively. Bounded motion occurs for 
Wf(£s) < E < Wj(^ u ); otherwise orbits are unbounded. Motion in the £,p$- 
plane is shown in Figure (6). The level curves of W% bear a strong resem- 
blance to the zero- velocity curves of the AKP in cylindrical coordinates (Fig. 
2). Whereas the latter curves represent motion boundaries for a given h z and 
E 7 the former correspond to actual trajectories. Lastly, we note that the equi- 
librium orbits derived in section 3, correspond to (£ S ,?7 S ) (stable, sombrero 
profile) and (£ U) ?7s) (unstable). 



5 Secular resonance 



The physical setting of the AKP is useful in identifying the dynamical re- 
sponse of a planet embedded in a accretion disk that sustains an asymmetric 
jet system. Such planets form inside the gas disk on circular planar orbits. 
As the disk's state of least energy is the sombrero profile (Paper II), the 
formed planet may assume one of the circular sombrero orbits and is there- 
fore immune to any additional dynamical excitation from the acceleration 
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of the star-disk system. In order to account for the large eccentricity of ex- 
trasolar planets through the AKP, the "initial" planetary orbital plane must 
be inclined with respect the direction of acceleration in order to display the 
librations described by the two Hamiltonians (Eqn. [T51 and ITT}) . The possi- 
ble physical realization of this geometry is when the jet system is inclined 
with respect to the star's rotation axis and therefore precesses around it. In 
this case, the acceleration has a constant vertical component that defines the 
sombrero profile and a precessing component that may excite the eccentricity 
of a planet embedded in the disk. This was the premise of Paper I where the 
excitation of eccentricity and inclination of extrasolar planets that formed in 
a jet-sustaining disk was examined. 

To illustrate the eccentricity excitation of a planet on an initially cir- 
cular planar orbit, it is instructive to first consider the situation where the 
precession frequency is so small that the acceleration appears to maintain a 
constant inclination with respect to the planetary orbit. In the context of the 
AKP, studying the excitation due to an inclined constant direction acceler- 
ation amounts to choosing initially inclined circular planetary orbits in the 
presence of a vertical acceleration. 

To examine this situation, the AKP was viewed in Paper I as a classical 
two-body Kepler problem perturbed by the acceleration A. Working in the 
secular regime where the excitation time \v\/A is much larger than the orbital 
period 27r/n, the acceleration potential R = A • x was averaged over the 
orbital period to find: 

3 3 

(R) = -tj ^- ' e = ~-^Aaesu\(m ~ J?)sini (21) 

where e and a are the orbital eccentricity and semi-major axis. The last ex- 
pression is obtained by choosing the z-axis along A; in this case, zu, f2 and 
I are the longitude of pericenter, the longitude of ascending node and the in- 
clination of the angular momentum vector h with respect to the acceleration 
A. The disappearance of the longitude in the perturbing potential implies 
that the semi-major axis is constant throughout the evolution. 

In terms of orbital elements, the conservation of h z can be written as 
Vl — e 2 cos I = cos Jo j where Jo is the initial inclination of the circular orbit. 
This allows us to reduce the problem to a single degree of freedom with the 
equations: 

. _ x/T - ? d(R) . _ VT~^ d (R) 



do; na 2 e de 



(22) 



where u> — zu — fi is the argument of pericenter and n = ^Jkj a? is the mean 
motion. The secular potential is now given as: 




(R) = --Aa\ - " esin^. (23) 



The orbital elements e and u> therefore follow curves of constant (R). In 
particular, the potential (R) (|23p has extremum points at: 



UJ 



±90°, e = v / 2sin(7 /2), and I = cos _1 ( v / cos7 ). (24) 
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This equilibrium may be interpreted as a secular resonance where the preces- 
sion rates of the pericenter and node become equal. The maximum value of 
e is sin Iq and corresponds to the cycle of initially circular orbits. The curves 
of constant (R) are shown in Figure (7) for an initial inclination of Iq = 45°. 

The value of (R) being constant throughout the evolution, it is possible 
to derive the analytical expression of the eccentricity evolution. Denoting by 
C the value of —2 (R) /3Aa, the equations ((22)) yield after some algebra: 



sin 2 Jo 



C 2 



c 1 



2 


\3A 


cos 


— t + cj) 




na 



+ 



sin 2 Iq 



C 2 



(25) 



where <f> is a phase corresponding to the initial value of e. The expression of 
sinw as a function of time can be deduced from the secular potential (|23p . 
For initially circular orbits, C = and u> = modulo 180°. In this case, 
equation l]25p reduces to the expression given by Equation (8) of Paper I: 



e(t) = 



3A 
2na 



sin ir 



(26) 



Equation (f2"5|) shows that the eccentricity time evolution occurs at the 
excitation frequency: 

n A = — , (27) 
na 

and that all libration cycles around exact secular resonance possess this same 
frequency. This was shown to be the case in Paper I where orbits were inte- 
grated numerically with the full equations of motion (i.e. without averaging) 
and a good match was found between the numerical evolution and equation 
(f2"6"| for initially circular orbits. 

To understand the origin of the secular resonance in terms of parabolic 
coordinates, we numerically integrated the equations of motion ([3]) for orbits 
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Fig. 8 Resonance orbits in parabolic coordinate phase space. The AKP parameters 
are those of a test particle orbiting a solar mass star accelerated by A = 5 x 1CP 11 
kms -2 . All orbits share the same initial osculating semi-major axis a = 1AU, 
argument of pericenter uj = 90°, longitude of ascending node Q = and true 
longitude 8 = 0. The different orbits correspond to different initial inclinat ions : 
Jo = 5°, 30°, 60° and 85°; their eccentricities are given by the expressions (I24| l. 
The equilibrium position along £-axis is larger for smaller To while the larger 7o, 
the larger the p,, ^-amplitude. The units are AU for £ and r/, and AUyear -1 for 
and p v . 



that were initially placed at exact resonance as defined by the osculating 
orbital elements given by The result is shown in Figure (8). It turns out 
that exact resonance orbits are on the stable equilibrium £ s in the £j>£-planc 
whereas the corresponding 7yp^-motion has a finite libration amplitude. This 
correspondence allows us to derive a geometric property of secular resonance 
orbits. As £ = 0, r + z (Eqn. [T^)) is constant for resonant orbits. This in 
turn shows that resonant motion is confined to the surface of the parabolic 
paraboloid: 



P 2 - n 2 
SO P_ 

2£o 



(28) 



where fo is the initial value of £. 
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6 Keplerian boundary of initially circular orbits 

The notion of a boundary in the AKP results from the possible equality 
at some relative distance of the two bodies of the gravitational interaction 
and the perturbing acceleration. Unlike the scale invariant classical two-body 
Kepler problem, the dimension of A naturally gives rise to a typical length 
that is proportional to (k/A) 1 / 2 . In Paper I, it was recognized that this length 
scale had to be an upper bound of the region of bounded motion and that a 
more precise estimate of the boundary is given by the equality of the mean 
motion n and the excitation frequency tia- This equality yields the limit 
beyond which the averaging procedure used in the derivation of (R) (Eqn. 
l2Tj) breaks down. The corresponding radius is given as: 





' k ' 


1 

2 


Okplr — 


M. 


~ 10 3 



2 x 1(T 12 kms" 2 " 


2 


' M ' 


A 




[m & \ 



AU. 



(29) 



In terms of ak p i r , the orbit that defines the outer edge of the sombrero profile 
(fTTTjl is given as: 



~crit 



akpir 



and p crit 



2V2 



akpir- 



(30) 



Using the formulation of the AKP in parabolic coordinates, it is possible to 
solve for the radius of the outer boundary of initially circular orbits whose 
plane normal is inclined to the direction of acceleration such as the orbits 
of the previous section and those considered in Paper I. These orbits are 
circular and planar but do not belong to the sombrero profile (i.e. are not 
least energy orbits) as they are initially inclined. For bounded motion to be 
possible, the Hamiltonian (|17p must have two equilibria. This in turn implies 
that (fT7|) has three roots at — for a given initially circular but inclined 
orbit. These roots are the solution of the equation: 



E = 



2e 



2 



(3-k 



(31) 



Calling a and Iq the initial osculating semi-major axis and inclination of the 
circular orbit and taking k = 1, the constants that appear in (|31[) are given 
as: 



1 

2a 



, h z = a cos Iq , A 



3o 



2 — , and j3 — — ' 

kplr 



6 V akpir 



(32) 



By substituting these constants into (f3"Tj) and writing a = (a/ak p i r ) 2 /6 and 
£ = S"ak P ir, the three roots may be obtained from the cubic equation: 



q £ 2 . COS 2 In 

aS 3 - — + (l- a)E — - = 0. 



(33) 
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The radius of the stable boundary is derived from (f33|) as the value of a 
for which the discriminant of the cubic equation vanishes; in this case, the 
equation does not admit three distinct real roots. This condition reads: 

a i , n a s a 2 9a sin 2 In sin 2 In 

a 4 - 3a 3 + (49 - 18 cos 2 I - 27 cos 4 J„) — + — — - - 0. 34 

v ' 16 8 16 

The solution of the quartic equation can be obtained by the usual standard 
methods; however, the expressions of the roots are quite cumbersome and are 
not given here explicitly. Instead, we use the solution of Eqn. (|34f directly 
to plot in Figure (9) the size of the stable boundary, a, normalized to akpin 
(or equivalently a) as a function of the initial inclination Iq. The region of 
bounded motion occurs inside 0.72 a^pir for all inclinations and its boundary 
decreases with increasing inclination, Jo- 
in Paper I, two examples of escape orbits were given: first, for ak p i r = 
100 AU an orbit was initially placed at an — 68.5 AU with an inclination of 
In = 30° (Figures 7 and 8 of Paper I). At such an inclination, Figure (9) 
yields a boundary radius of an = 68.5 AU. The chosen orbit was therefore 
on the unstable equilibrium point of the £-Hamiltonian. The second exam- 
ple was that of the stellar binary companion to v Andromedae. To illustrate 
how to reach the current orbital configuration of this multiplanet binary sys- 
tem, a time dependent acceleration with a minimum a^pir = 300 AU was 
applied to the configuration where the stellar companion was initially placed 
at a = 298 AU with In = 60°. In this case, Figure (9) shows that the min- 
imum stable boundary radius corresponds to a = 198 AU. As the accelera- 
tion pulse was reaching its maximum, the stable boundary radius decreased 
from infinity and approached 198 AU. The stellar companion was therefore no 
longer gravitationally bound to v Andromedae and moved away from it (Fig- 
ure 14 of Paper I). However as the orbital time of the companion's motion, 
~ 5000 years, was larger than the acceleration decay timescale of 2000 years, 
the stellar companion was not left with enough time to escape before the 
boundary radius bounced back to infinity. The net result of the accelera- 
tion led to the radial migration as well as the eccentricity and inclination 
excitation of the stellar companion. 



7 Further developments 

In this paper, we have examined the dynamics of the accelerated Kepler 
problem in detail using the problem's integrability. We have (i) compared the 
AKP to the classical Kepler problem, (ii) derived analytically the features of 
least energy orbits, (iii) derived the general solution of the secular problem 
and interpreted the secular resonance in terms of parabolic variables. Lastly 
(iv), we used the formulation of parabolic variables to determine analytically 
the boundary of initially inclined circular orbits that are of particular impor- 
tance to the problem of radial migration in binary systems as well as to the 
truncation of accretion disks through stellar jet acceleration. 

The AKP's integrability makes it an interesting Hamiltonian system that 
can be perturbed by additional processes. The perturbation may result in the 
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Fig. 9 The size of the Keplerian boundary of initially circular and inclined orbits 
in the AKP. The ratio of a to a^ p i T is shown as a function of Iq (degrees). 



loss of integrability and may lead to possible diffusion. In turn, such mech- 
anisms may further enhance eccentricity excitation and radial migration. In 
the context of the dynamics of forming protoplanets in a jet-sustaining disk, 
the modeling of a general time dependence of the perturbing acceleration 
provides a more realistic model of jet-disk interaction with further conse- 
quences such as disk heating (simples examples are given in Papers I and 
II). 

An interesting simple dynamical problem that currently remains unsolved 
is the possible integrability or quasi-integrability of the precessing AKP. The 
resonance studied in section 5 illustrated how to excite the eccentricity of an 
initially circular orbit in the case where the constant-amplitude acceleration 
has a small precession frequency compared to the excitation frequency ua- 
The more general case of a precessing constant-amplitude acceleration (that 
corresponds to a precessing jet) yields the precessing AKP equations: 



dv 

dt 



du 

Au, and — = it a x u, 
dt 



(35) 



where £Ia denotes the acceleration's constant rotation vector, A is constant 
and u is the precessing direction of acceleration. To recover the setting of the 
AKP, two rotations may be applied to ([33)1 . First, a rotation at the rate Qa 
around £1 a rewrites the equations in the rotating frame where the z— axis is 
chosen along 17^. In this rotating frame, the acceleration (or the vector u) 
is confined to the ccz-plane and makes an angle a with the z-axis. Second, 
a rotation of angle a around the y-axis makes the direction of acceleration 
along z. The corresponding transformations are simple to carry out and the 
resulting equations are: 



X — 2Qa cos a Y = fi A (cos 2 a X + sin a cos a Z) 

kX 



(X 2 + Y 2 + Z 2 ) 3 / 2 



Y + 2f2 A (cosa X + sina Z) = Q A Y 



kY 



(X 2 + Y 2 + Z 2 fl 2 ' 



(36) 
(37) 
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Z — 2f2 A sin a Y = fi\ (sin a cos a X + sin 2 a Z) 

- (x 2 + r 2 + z 2 ) 3 / 2 + A (38) 

With respect to the AKP, precession introduces the Coriolis force (second 
term on the left hand side of each equation) and the centrifugal force with 
respect to the rotation vector Sl A (first term on the right hand side of each 
equation). The momenta of the precessing AKP are given as: 

p x = X - fl A Y cos a, (39) 
p Y = Y + Q A (X cosa + Z sin a), (40) 
pz = Z- Q A Y sin a, (41) 

and the corresponding Hamiltonian is: 
p 2 

H = — + Q A Y(px cos a + pz sin a) 

k 

- Q A p Y {X cos a + Z sin a) — — — — AZ. (42) 

l x l 

This in turn shows that the energy integral of the AKP is: 

E = +r +Z - tfA [(Xcosa + Zsina) 2 + Y 2 } - ±- - AZ. (43) 

As axisymmetry has now been lost, the vertical component of angular mo- 
mentum is not conserved. We have searched for analytic integrals of the 
precessing AKP without much success. We note however that averaged pre- 
cessing AKP is integrable and that an exact analytical solution was derived 
in Paper I. Further numerical investigations will shed some light on the in- 
tegrability of the precessing AKP. 
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